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There exist a wide temperature region (GiT < T — T c < TyGi), where the 
influence of fluctuations on the thermodynamic properties of superconductors can 
be taken into account in the linear (Gaussian) approximation, while their influence 
on the kinetic properties is strongly nonlinear. Maki-Thompson cotribution to 
conductivity saturates in this region. However, Aslamazov-Larkin contribution 
becomes even more singular. Such an enhancement is connected with the fact that 
nonlinear effects yield an increase in the lifetime of fluctuating pairs. Pairbreaking 
and energy relaxation processes can diminish the nonlinear effects. 

1. Introduction 

Electron scattering off the usual impurities leads to a temperature inde- 
pendent residual resistance of the normal metal [1]. The conductivity of 
bulk samples and films can be measured with a very high accuracy. This 
gives a possibility to study different mechanisms, leading to the temperature 
dependent conductivity at low temperatures. One of these mechanisms is 
connected with the thermal fluctuations above the superconducting transi- 
tion temperature T c [2-5]. There are two kinds of fluctuation corrections, 
leading to the temperature dependent conductivity above T c . The first one 
is known as Maki-Thompson contributions and the second one is conduc- 
tivity of fluctuating pairs (Aslamazov-Larkin contribution). These correc- 
tions have different dependence on the spin flip scattering time t s . Charac- 
teristic temperature range for contributions of both types is Ginzburg pa- 
rameter Gi, where Gi depends on dimensionality and for films is equal to 



Gi = ro = l/32uDd = e 2 /16ha n - Where v = mp 2 /2ir 2 is the electron den- 
sity of states per spin, D = \vphr is the diffusion coefficient, d is the film 
thickness, l tr is the electron mean free path, p is the Fermi momentum, cr n is 
the conductance of a square film. It has been found in paper [6] , that nonlin- 
ear fluctuation phenomena lead to a new temperature scale T c \/Gi (see also 
[7-10]). In this paper we obtain expressions for conductivity in the temper- 
ature region Gi < r < y/Gi, where both the Gaussian approximation works 
well and the nonlinear fluctuation effects are important. 

In paper [6] an attempt to find the fluctuating correction to conductivity 
was made. The main point was that long wave fluctuations with Dk 2 < Tt 
are essential. Such fluctuations can be considered as a Bose condensate. 
Dynamics of superconductors should be considered on the background of 
these fluctuations. They give a pseudogap in the excitation spectrum. In 
this paper we will show, that short wave fluctuations with Dk 2 ^> Tt can 
be important. It was found in paper [11] that short wave fluctuations of 
the order parameter A act on the electron Green functions as paramagnetic 
impurities with depairing factor Y = t s -1 = (\A\ 2 )/e. Essential values of 
energy e are of the order of e ~ A ~ Ty/Gi. Hence V is of the order of 
Ty/Gi. So large value of the depairing factor leads to saturation of the M.-T. 
contribution to conductivity in the temperature region r < v^Gi. 

A more complicated situation takes place for A.-L. contribution. This 
contribution is proportional to the density of pairs and their lifetime. For 
large enough values of r the time dependent Ginzburg-Landau equation 
(TDGL) can be used to obtain this lifetime. It is proportional to h/(T — T c ) 
and hence A.L. contribution proportional to r _1 . If the concentration of 
paramagnetic impurities is large or if the energy relaxation time is short, 
the TDGL equation can be used for all temperatures T. In this case A.L. 
contribution is valid in temperature range r > Gi. But in the opposite 
limiting case, the nonlinear fluctuating effects destroy the applicability of 
TDGL equation and lead to the increase of lifetime of fluctuating pairs. As 
a result the A.L. contribution to conductivity becomes more singular in the 
temperature region \/Gi > r > Gi. 
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Qualitative picture 



In the temperature region 1 3> r ^> Gi thermodynamic fluctuations of the 
order parameter A can be considered to be Gaussian. The corresponding 
correlator has the form: 

/a*a - T 1 _ 256 GiT 2 

In order to calculate thermodynamic quantities in the temperature region 
t > Gi it is sufficient to know this correlator only. However, to calculate ki- 
netic coefficients a more complicated problem has to be solved. One has to 
find how Gaussian fluctuations change the one-particle excitation spectrum. 
The longwavelength fluctuations with k 2 < k^ a = 8Tt/ttD can be consid- 
ered as a local condensate. They lead to the formation of the pseudogap 
in the one-particle spectrum of excitations. From Eq.([TJ) we see that the 
pseudogap is equal to 

A PG = -Voir. (2) 

7T 

Not very close to the transition (r > v^Gi) only the excitations with energy 
uj > Ape are important. The pseudogap does not play any role for such 
excitations. Thus, it is sufficient to consider fluctuations in the linear ap- 
proximation only (see [3], [4] and [5]). It is important, however, that in the 
temperature region r < \^Gi the excitations with energy to < Ape become 
essential. In [6] the fluctuation correction to conductivity was considered in 
view of the pseudogap. The pseudogap was taken into account in the same 
way as the gap below the transition temperature. This approximation gives 
a correct estimate for the width of the temperature region where the non- 
linear effects are important. However, the model considered in [6] with A 
being a constant can not reproduce the correct temperature dependence for 
conductivity in the temperature region r < \^Gi. 

To describe the nonlinear effects we will consider fluctuations of A in 
the statical approximation. This is eligible, since the lifetime of fluctuations 
(Tr) _1 is large compared to the inverse pseudogap. However, the spatial 
dispersion of the pseudogap changes the physical picture significantly. To 
take into account the spatial variations, we have to calculate the conductivity 
as a function of the order parameter A(r) being an arbitrary function of r and 
average out the result over the Gaussian fluctuations with the correlator ([!]). 
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In the present paper we accomplish this program up to a numerical coefficient 
in the limiting case when the energy relaxation rate is large (r £ <C (Tr) -1 ). In 
the other cases, we obtain a functional form of the temperature dependence 
of conductivity without calculating coefficients. 

To consider the spatial dependence of the order parameter we will use the 
results obtained in [11]. In this work, it was shown that the spatial variations 
of A act on one-particle excitations in the same way as do the magnetic 
impurities. In this case, the total pairbreaking rate T can be written as a sum 
of the pairbreaking rate due to the magnetic impurities and the fluctuation 
term. Thus, the self-consistent equation for T takes the following form 

r d 2 k (A*A fc ) 1 

J (2vr) 2 u + \DW + T^ r s - 1 ' 

It is important to mention, that Eq.([3D is exact either if u 3> T or if r s is 
very small, so that the first term in Eq. (|3|) is a small correction to the second 
one. In the other cases the self-consistent Eq.(|]) can be considered as an 
estimate and gives the result valid on the order of magnitude only. 
In the region u < T, T > Tr we obtain from Eqs.(|3]), (|1|): 

r^fon^f , (4, 



7T 



which coincides with the value obtained in [7,12] up to the logarithmic term. 
Below, we repeat the derivation done in [11] and show that the pseudogap 
does not change result (Q) qualitatively. 

Let us note, that the pair-breaking rate Y is of the order of the pseudogap 
Ape- Thus, a wide maximum appears in the density of states at to ~ Ape- 

As it is known from [5], the Maki-Thomson correction to conductivity 
saturates for Tr < T and takes the form: 

6a MT 8TGi vrr 

= — F^ ln 77^- (5) 

<Jq 7rl 41 T 

As it can be seen from Eqs.(^J^) such a saturation takes place when r < 
v^Gi. The similar results have been obtained in [7,8,10]. However, numerical 
coefficients are different. 

Note, that the numerical coefficient in Eq. (H) depends on the way how the 
summation of higher order diagrams is made. However, its exact value is not 
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very important since in the region Tr < T Maki-Thomson contribution is less 
singular compared to Aslamazov-Larkin contribution and can be neglected. 
Let us realize, that AL contribution does not saturate when T tends to T c 
but becomes more and more singular. 

To estimate AL contribution due to the appearance of the fluctuating 
Cooper pairs one can use simple Drude formula: 



ne 
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5a AL = ^-r f , (6) 



in 



Where n, m and Tfi are concentration, mass and lifetime of the fluctuating 
Cooper pairs. The ratio n/m can be estimated from Eq.(|l|), while the lifetime 
follows from the TDGL equations. 



d_ 



(^ + rf + -Tr A 4 (i)=((t), (7) 



where ( is the Langevin noise. In two-dimensional case we have 

n T 



m 2ndh 

and 



8(T - T c )a 

Not very close to the transition (Tr > Ape) or if the energy relaxation rate 
is very large we can put a — 1, since the quasiparticles are in the thermal 
equilibrium. Thus, we have 

6*AL = Gi (g) 
a t 

In the presence of the pseudogap there is no equilibrium and the coefficient 
a becomes greater then one. Recall, that analogous changes in the coefficient 
a in TDGL equations for |A| appear below the transition temperature (see 
e.g. [13,14,15,16,17]). The growth of the coefficient a and, consequently, the 
growth of the fluctuations lifetime is due to that the quasiparticles require 
more time to attain the thermal equilibrium (we denote the corresponding 
time as r e ). A rough estimate gives a ~ ApcT e . In the case of weak energy 
relaxation, r e has to be determined from the diffusion equation in view of 
pseudogap (see [18,19,20]). Note, that in this complicated case coefficient a 
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becomes a non-local operator. Rough estimates give the following value for 
the thermal equilibrium transition time r e ~ (Dk 2 ^^)^ 1 ~ (Tr) _1 . Taking 
into account Eq.(|2|) we obtain 

^ Gi 3 / 2 

— = —T- 9 

(Jo ' 

We see, that paraconductivity can exceed normal conductivity <To in the 
region Gi 3 ^ 4 > t > Gi. Let us emphasize, that in this region corrections 
to all the thermodynamic coefficients are small and the linear theory is well 
applicable for them. 

Let us now discuss the role of the energy relaxation processes, character- 
ized by a quasiparticle lifetime t £ . Non-elastic electron-electron scattering 
in dirty metals leads to the following electron-electron collision times in 2D 
case 

t~ x -Tdlp 2 ~ GiT. 

Such a large collision time does not change nonlinear effects. However, 
nonelastic electron scattering off phonons and other possible collective ex- 
citations can decrease r e significantly. These processes together with addi- 
tional pairbreaking processes (due to magnetic impurities or magnetic field) 
lead to a decrease of the nonlinear effects. Energy relaxation reduces the 
thermal equilibrium transition time r e . If these processes are very strong 
(for example if the temperature is relatively large) transport equation for the 
distribution function becomes local and in the limit Tr ~ Dk 2 <C r" 1 we 
can write r e = r e . Thus, in the temperature region under consideration we 
have 

8a G^Tr r 



1 e 

C T 



Elastic scattering off magnetic impurities and magnetic field also tend to di- 
minish the nonlinear fluctuation effects in conductivity but in a different way. 
These scattering processes (as well as scatterings off the static fluctuations 
of the order parameter) do not affect quasiparticle motion and r e thereby. 
However, if the pairbreaking rate is large enough (r > A PG ) these processes 
lead to the reduced pseudogap A PG ~ (|A| 2 )/T (recall, that without pair- 
breaking we would have Ape ~ \/(|A| 2 ) ~ TGi 1/2 ). Thus, the fluctuation 
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correction can be written as 



In the presence of both a strong pairbreaking and a large energy relaxation 
it is possible to derive exact expressions with logarithmic accuracy for the 
coefficient a in TDGL, which is local in this case, and for paraconductivity. In 
this case, the main contribution to a comes from the fluctuations with Tr < 
Dk 2 < TgT 1 . The first inequality allows us to consider only the leading terms 
in the expansion a over A, while the second implies a local approximation in 
the transport equation. The result is 

« = ^, (12) 



5a 32Gi 2 T 2 r ( 



7T 



In . (13) 

a 7r 2 IY 8Tt £ t y ' 

Note, that Eqs.(|9|-|r3"D are valid only if the parameters T and r e are such that 
the contribution to conductivity 5a is larger than the usual Aslamazov-Larkin 
contribution Eq.(H). If T > T, Tt £ < v^Gi or if T 2 r £ /T < Gi, than nonlinear 
effects are negligible and the usual result is valid for all r > Gi. Note, that 
MT contribution saturates at temperatures when Tr ~ max 



i TVGi 



2. Depairing factor induced by fluctuations 

Nonzero fluctuating order parameter A and Gor'kov green function (3 [6] exist 
above the transition temperature. In the temperature region r > Gi the main 
contribution to the value of order parameter A arises from zero "frequency". 
The momentum space can be separated into two parts: 7rDk 2 /8T < r and 
7iDk 2 /8T > t. The fluctuations with irDk 2 /8T > r can be considered 
as "fast" variables, created on the background of slow fluctuations with 
7iDk 2 /8T < t. The "fast" fluctuations induce intrinsic depairing factor T, 
even if external depairing factor, connected with paramagnetic impurities is 
missing (r s — ► oo). Similar phenomenon was studied in paper [11]. Using 
the method, developed in this paper, we obtain expressions for the statical 
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Green functions a, (3 and depairing factor T. We start from Usadel equation 
for Green functions a and j3 in the dirty limit (see [6,21]): 



D 



Aa-uf3+ — (aV 2 (3 - (3V 2 aj = a/3T. 



(14) 



According to paper [11] we present Green functions a, (3 in the field of 
"fast" fluctuations of the order parameter A(k) in the form 

a=(a) + a 1 , P = (P)+Px (15) 

The deviations of the Green functions from their mean values can be 
found on the perturbation theory. We have [11] 

A(k)(a)((3) 



(16) 



u(a) + (A)((3)+Dk 2 /2 

The "mean" Green functions (a), (/3) are solutions of the following system 
of equations: 



(a)(A)-u(P) = (a)(P)T 



(«) 2 + (P) 2 = ] 

The value of parameter T is defined by Eq.fllDD and equals to 

d 2 k (A* k A k ) 



where (A) 



(2tt) 2 (a)w+(A)(P) + Dk 2 /2' 
(\A\ 2 y/ 2 . The quantity (A) in Eqs 



(17) 



(18) 



E3) should be 

understood as an integral over k from the expression (|18D in the range 
nDk 2 /8T < t and takes the form 



(A) 



T 
ud 



d 2 k 



(27r) 2 r + f^A; 2 



,1/2 


r64Gil 


« T 




7T 2 



From Eqs(|T|) and (|18|) we obtain 
„ 16TGi 1 



In 



TIT 



' n(co(a) + (A)((3)y 
ATt 



(19) 



(20) 



£-(u(a) + (A){P))-l 

As can be seen from Eq. (p0|) , T(u) is a function of energy uo. In the range 
r < -\/Gi essential values of u are of the order of T. Thus, T itself is of the 
order of (A) (see (19)). This order of T is connected with fluctuations of the 
modulus of the order parameter. This value is much larger, then the one due 
to the phase fluctuations of the order parameter (see [6]). 
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3. Equations for the time dependent order pa- 
rameter 

Statical Ginzburg-Landau equations are valid in wide temperature region 

Gi« \1-T/T C \ < 1 (21) 

TDGL equations are valid if the energy relaxation time r £ or the pair 
breaking time r s = T _1 is short enough [13-16]. For large r e 's the dynamics 
of normal excitations becomes essential. As a result, the dynamical term in 
the equation for the order parameter becomes more complicated. Below, we 
derive the corresponding equation. 

Order parameters A 1)2 (i) can be written as 

Ai, 2 (0 = ^if 2 ( M ) (22) 
where Green function G can be presented in the form [18] 

6 =(T'T*) (23) 

where G R,A ' K are retarded, advanced and Keldish Green functions. Each of 
them is matrix of Gor'kov-Nambu 

G R,A,K_(gi] A-( ° ] Al \ (24) 

G - { -F 2 ; g 2 ) ' A -{ -A 2 ; ) ^ 

where A 2 (cu) = A* (—a;). 

In the dirty limit we have the following system of equations for G R ' A (see 
[19]) 

- Dd T (g R > A d^ A - F.Y^pj + 2iA h2 g R ' A - (25) 
-2ieF^ A + -g RA F^ A = _/f 

T s 

where the electron-phonon collision integral j P ^ R ' A ^ j n the vicinity of the 
transition temperature T c for small values of energy \e\ <T is equal to 

jPHR,A) = ± L Fi y (26) 



The Keldysh Green function G K can be presented in the form [20] 

G = J dhiGXf-fG*) (27) 

where distribution function / is given as [20] 

/ = / + r z h (28) 

Equations for the distribution functions has been derived in paper 
[20] and have the form 

-4r{%^- GRG ^}- D l^ + 2 % Spa+ < 29 > 
+4/ 2 p "(/,)=0. 

where 

j e = SprJG R dG R - G A dG A ), d=^--ieAr z , (30) 

or 

2a = G r t z - t z G a , 25 = G R - G A , 2 7 = G R + G A 



In the important limiting case £~T> A the Eqs. (|25|j29|) can be simpli- 
fied and we obtain 



^ = - i [r T fe-f| 5 ] A; (31) 
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+^ i (F 1 R -^)} + /f"(/)=0 

f Wt + 1 ("^ + ^) + ^ + *>) } + = 0, 

where 

j E = -F*^l + F 2 R ^l + Ff^l - F A ^- (32) 
dr dr dr dr 

The general expression for collision integrals I±% is given in papers [20,22]. 
For small values of energy |e| <C T these integrals can be taken in the simple 
form 

'**</>=£(-*(£)+/) (33) 

r- 1 = 7C(3)7r^ 2 T 3 /2( S p) 2 

where s is the velocity of sound in metal and g is the electron-phonon coupling 
constant. 

In the limiting case of strong energy relaxation, when r £ A <C 1 the dis- 
tribution function / can be taken as the equilibrium one 

/ = th(e/2T), h = (34) 



In this case one can obtain from Eqs.(^), ( |3~TD and (33) time-dependent 
Ginzburg-Landau equation in the usual form 

( 1 - T/T c - | A\A A + ^ A - ( £ + 2ie^ A = (35) 

\ / 87r 2 T 2 ' 1 ; 8T - 8T v<9t v v ; 

If condition r £ A <C 1 is not satisfied, then the deviation of the distribution 
function / from it's equilibrium value can lead to the change of last term in 
Eq.@. 

In the range T ^> A the crossing term in Eq. (|3T|) is small on the parameter 
(A/r) 2 . Thus, in the leading approximation system (31) is diagonal. 
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With the aid of Eqs.fl27|), (f28|) and (§1]) we can rewrite Eq.(|22|) in the 
form 




-iui — D 



Q r 2 



T + 



7T 

8T 



— iuji — D 



dr 2 



In Eqs.(|36|) we put 



A r 



./ 



7T 

"2 

7T 
2 



de 
2tt 

cfe 
2^ 



8f{F*-F?)-f x (F« + Fi 



R 



6f(F 2 R - F 2 A ) - h 



= 

(36) 
= 

(37) 



th(e/2T) + 5f 

In ( ]36| ) the contributions due to the second terms are of the order (A/T) 2 . 
This result is due to the cancellation of terms, going from quantities 5f and 
j\. But in the next orders of the perturbation theory the quantity f\ becomes 
small and the main contribution arises from distribution function Sf beyond 
the perturbation theory. 



4. The conductivity of fluctuating pairs 
(Aslamazov-Larkin contribution) 

The conductivity of fluctuating pairs is given by diagrams on Fig. la. We 
suppose below, that order parameters Ai^ can be written as a sum of two 
terms. One of them is connected with the statical thermodynamic fluctua- 
tions A, A*. In the range r > Gi these fluctuations are Gaussian with the 
correlator given by Eq. ([!]). The wavy line on Fig. la gives the dynamical 
fluctuations A 12 of the order parameter. The correlators of these fluctuations 
Kij should be found on the background of thermodynamic fluctuations 

Kvfa) = uiAtAjU (38) 

The contribution to conductivity can be expressed through correlators K 
in the same way as in the case of weak fluctuations [3]. 

First, we have to find conductivity as a function of Matzubara frequency 
u and, then, perform analytical continuation on ujq. The correction to cur- 
rent was found in paper [6] with the aid of the equations for Green function 
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Fig. la,b a - Aslamasov-Larkin contribution to the conductivity 
b - Maki-Thompson contributions to the conductivity 



in the dirty limit in high frequency fields 

J: = Y d j d^T^SpL^K^ + uj^r.n^K^n.r^ (39) 



where A m is the vector potential of the external field, and matrix L is equal 
to 

ireD d - 



-^12 — -^21 — 0) -^11 ( r ) 



t a r a 

-^22 — -^11 



2T dr c 

After the analytical continuation over uj in Eq.(R3) we obtain 



(40) 



1 

2d 



1 



1 



uo\ — iu — 5 + 5 



s P [L a r k( 



UJ\ — IUJ 



+ S,r,r 1 )L^k(u l -8, 



ri,r 



At 



(41) 



It was found in paper [6], that fluctuations are weak in the range r > 
(Gi) 1//2 . And in this region we have 



KniuJt + 5) = K 22 (uj 1 + 6) 



(42) 



From Eqs. ([E|-|42|) we obtain the well known result for the paraconductivity 
[3] 

a^/a = ^ (43) 
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To obtain the conductivity in the temperature region r < (Gi) 1 / 2 we 
should find correlation functions K in the field of thermodynamic fluctuations 
A. Than we have to average out over A the expression for conductivity The 
correlation functions K can be found from Eq.(p6|). 



K' 



T + 



sr 
C21] 



c 



11, 



-c 



12 



T 



where operators Q,- are given by 



C 



11 



7T 

2 



27 



Ff-FtfSfW-fFlFf + Fi 



(1) _ f (!) 



c. 



22 



(44) 



(45) 



C 



12 — 7T 



21 



22 



7T 

2 

7T 

2 



-00 
00 



de 
27 

27 

de 
27 



F A 



F 



F 1 A )6fW-fr(F? + Fi 



(2) 



5fW + (F 2 R + F 2 
SfM + (F« + F 



(i) 



(2) 



In equations (45) the operators 5f^ 1,2 \ /1 are defined as 

5f = 5f^A 1 + 5f^A 2 , (46) 

/1 = A (1) A! + A (2) A 2 

and 5/, fi are solutions of the system (31) in field of Ai, A 2 . For an arbitrary 
function A(r) system (31) can not be solved analytically. Nevertheless, in 
the range r < (Gi) 1 / 2 the expression for correlation functions K can be found 
with the logarithmic accuracy if the value of the external depairing factor V 
is larger then A. In this case simple expression for Green functions F^ 2 ' can 
be used 



R,A 



-iA 



iR,A 



-iA* 



(47) 



r^ie' 2 T^ie 
If Dk 2 |A| 2 /r, the contribution of Sf^ 1 ' 2 ^ is cancelled out in quantities 
C u and C 22 . Note, that if Dk 2 < |A| 2 /r, than f[ h2) < 5f^ 2 \ Thus, this 
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region gives the dominant contribution to C^. Eqs.([4"4",l"5) in this case can 
be reduced to the following form 



8T 



K ll + ^^ + r->-D^y\A-K u ) + (48) 



7ru;iA 
167T 



TTUJiA 

+ 16TT 

t+— (uJt-D— 
or z 



UJl + T e 
2 



K 



12 



162TV ^ e dr2 

-D*- 

dr 2 

f 16TT 



) 1 (AK 21 ) = 5(r-r 1 ), 



d 2 N -i 



( Wl + r- 1 - D 



d 2 
dr* 



D W 2 > i A *^ 12) + 



22 



0. 



8T 



i^2i + 



i6Tr 



(o; 1 + r e - 1 -^)" 1 (A^ 21 ) + 

<9 2 x- i 



I — (wi + -r 

1 RTF \ 1 1 



16TT 

T+—(UJi-D — 



-1 



T> 



A*Ky 



0. 



8T 



<9r 2 



^22 + 



vr^iA* 

i6Tr 



-i 

£ 



n^A* 



a 2 s-1 



16TT 1 ,£ Or 2 ' 
This system can be solved with logarithmic accuracy in the case of strong 
energy relaxation r^ 1 > Tr. In this region we obtain from Eq.(|T]) 

d 2 



A) 



64Gi 



TV In 



7T 



7T 



(49) 



From Eqs. ([I]), fl4"8] ) and (^5) we obtain the following expression for the 
correlators K: 

t + —k 2 H ^ In ( 



8T 

2 / 4GiTr £ o; 1 " 



7lT 



I\K 



ii 



1Ttt f 



K 22 = K u 



(50) 



where 



dxdy 



(x + l)(y+ l)yj(x - y) 2 + 2(x + yja + a 2 



(51) 
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1 , ^Tr.Gi , / vr \ 
o = H In — 

tiTt \8Ttt £ J 

The nondiagonal elements in K give logarithmically small contribution 
to conductivity. As a result we obtain 



a a 32Gi 2 T 2 7> 



0o 



7T 



In 



7T 



(52) 



The situation becomes more complicated if energy relaxation time r £ is 
large. From ( [48]) we obtain the following equation for correlator K\\ 



7T / or 



d 2 



TlLJiA I i 



(53) 



7ro;i 
16TT 



i6rr 



<9 2 



D- 



D 



A 



dr 2 

-i 



A 



7T 

8T 



A* (oji + t~ 1 — D 



LVi-D 

d 2 



dr 2 

-i 



+ 



A*K U 



= 5{r — ri) 

First, we find the mean value for the following quantity 

2 



(A'fa-D— y\Aexp(tkr))) 



(54) 



64T 2 Gi 



' vr (Dfc 2 + ^) 2 ^ 



2 Dk 2 + ^ 



In I 
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From here we see that the coefficient near U\ in equation for Ku is log- 
arithmically large. Contrary to the previous case (r^ 1 ^> Tr) the last term 
in the right-hand side of Eq. (|53|) is essential and leads together with nondi- 
agonal elements in K to the cancellation of large terms in conductivity. To 
check this, we have to find the mean value of the product of four A in last 
term in Eq. (|53D . We have 



h 



7T 



i6Tr 



A 



7T 



d 2 



T + it* K - + ( 55 ) 



8T 



dr 2 
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TTUJl 
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exp(i/cr) 



A* [ui - D 
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-1 
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-i 



A* [u! - D 
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y r d%_ r d 2 k 2 
lQTTud j J (2tt) 2 J J2t7) 2 



(wi + D(k - k 2 y){LOi + D(k - hY) ( r + —Dk 2 3 + u) X a k 



ST 
dr 2 

irD 

' w 

7i 



A* exp(ifcr)) 



-i 



-i 



where 

k 2> = k — k 1 - k 2 , 



4TGi 



Oik 



In 



vrr Dk 2 + ^ V 

7T \ 1 



(56) 



The In 2 term can be easy separated from expression fl55l) . As result we 
obtain 



r + ^fc 2 + o;ia fc 



47r«fcGi 



7rD 

g(fc 2 - A; 2 ) +cui(a fcl - gfc) 



(57) 



(2vr) 2 ( T + m k \) (A* - kf (r + - k) 2 + uj x a k ^ k ) J 

In Eg . (|55"D we omitted the "diagonal" term with denominator of type 
[u\ + £)(/c + k{) 2 ] 2 . This term leads to a small correction to coefficient near 
uj\ in (|53|). 

With the same accuracy we present here the expression for nondiagonal 
elements Ki2,K 2 \ as 



K> 
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l&TT 



d 2 k\d 2 k 2 



(58) 



A* kl A* k2 Kn(k) 



{uJi + D(k - h) 2 ) t + ?§{k -h- k 2 ) 2 + WKXk-k^ 



K 
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16TT 



d 2 k 3 d 2 k 4 



A fc3 A fc4 X 22 (fc) 



(wi + D(fc + h) 2 ) t + mk + h + k 4 ) 2 + Lu iak+k3+k4 
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With the aid of Eqs.( |5T| , |5"8|) we obtain that correction to conductivity in 
the form 



a a 4TGi 2 



(59) 



do ttTt 2 

This expression is valid up to a numerical factor of the order of unity. 

If the external depairing factor T is zero (superconductor without para- 
magnetic impurities), then the quantity T in Eqs. (|51|j59|) should be replaced 
by its intrinsic value 

T « TGi 1/2 (60) 

(see Eq.(0)). 

As a result we obtain in the temperature region Gi < r < (Gi) 1 / 2 

- - ( 6l ) 

cr TIT 

Eq.(^) means, that A.L. contribution to conductivity is strongly en- 
hanced in the temperature region Gi < r < (Gi) 1 / 2 . 

5. Maki- Thompson contribution to conductiv- 
ity in nonlinear on fluctuations region 

The general expression for Maki-Thompson contribution to conductivity (a b ) 
was given in paper [6]. Equation (28) in [6] can be considered as interpolation 
expression for Maki-Thompson contribution to the conductivity, that is valid 
in the entire temperature region r > Gi. The depairing factor T in Eq.(28) 
in paper [6] should be changed to be a sum of two terms: the first one is 
the external depairing factor t~ x , connected with the spin flip scattering on 
magnetic impurities, and the second one is the intrinsic depairing factor, 
given by Eq. (|20"D . As a result we obtain 

a b 7T r d 2 k 1 1 

V ~8d^J (2vr) 2 r + Dk 2 /2 TTfp 2 ~ (M 

2Gi 



T 
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In the range Gi < r < (Gi) 1 / 2 Maki-Thompson contribution reaches its 
saturation value and effectively becomes temperature independent 

^ = (Gi)1 / 2ln (M^ (63) 

The correction remains small in the entire region where nonlinear effects 
are important Gi < r < (Gi) 1 / 2 . 

Note, that real superconductors are always inhomogeneous. The finite 
value of the transition width leads to the appearance of an effective depairing 
factor [11]. And the value of this depairing factor can be large enough in the 
units of TGi. In such a case Maki-Thompson contribution to conductivity 
is small compared to Aslamazov-Larkin contribution in the full temperature 
region. 



6. Conclusion 

We see, that nonlinear fluctuation effects in kinetics phenomena are much 
stronger than in thermodynamics. If external depairing factor is absent, 
then the nonlinear effects lead to the saturation of Maki-Thompson con- 
tribution to conductivity in the temperature region r < (Gi) 1/2 . In this 
temperature region Aslamazov-Larkin contribution becomes even stronger 
and grows as <J a /a ~ Gi 3 ^ 2 /r 2 . In a superconductor with a large enough 
external depairing factor T = t s _1 > T(Gi) 1 / 2 or short energy relaxation time 
tT 1 > T(Gi) 1 / 2 the Maki-Thompson contribution saturates in temperature 
region Tr < T or Tr < r £ _1 . The nonlinear effects are not very important 
for it. Magnetic impurities and energy relaxation act on Aslamazov-Larkin 
contribution in a different way. Energy relaxation leads to the appearance 
of collision integral in the kinetic equation for the distribution functions of 
normal excitations. This collision integral diminishes nonequilibrium con- 
tributions to the distribution functions. Magnetic impurities and magnetic 
field act only on the superconductivity and don't lead to the relaxation of the 
distribution functions. However, TDGL equation depends essentially on the 
electron distribution function. If parameter t" 1 > T(Gi) 1 / 2 , then nonlinear 
fluctuation effects are not essential and Aslamazov-Larkin contribution re- 
mains the same o a /<j = Gi/r in full temperature region r > Gi. If inequality 
t~ 1 < T(Gi) 1/2 is fulfilled, then the law a a /a ~ [Gif^jr 2 takes place in 
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the temperature region Tr > t~ x . In the region (Tr £ ) _1 > r > Gi correction 
to conductivity is given by the expression a a /cr ~ Gi 3 ^ 2 Tr e /r (see Eq.(|52])). 
Magnetic impurities (or current) suppress nonlinear fluctuation effects in a a , 
but the effect is not as strong as in the case of energy relaxation. In the 
range TGi/T > r > Gi the correction to conductivity a a is given by Eq . (|59|) 
(T a /a ~ TGi 2 / (Tr 2 ). In the temperature region r > TGi/T the correction 
a a has the form (^) in the linear approximation. 

It is essential, that conductivity of fluctuating pairs can be larger then 
conductivity of normal electrons in the temperature region, where correction 
to the thermodynamic quantities is still small (see Eq.(61)). 
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